Intermittency route to chaos for the nuclear bilhard - a quaUtative study 
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We analyze on a simple classical billiard system the onset of chaotical behaviour in different 
dynamical states. A classical version of the "nuclear billiard" with a 2D deep Woods-Saxon potential 
is used. We take into account the coupling between the single-particle and the collective degrees 
of freedom in the presence of dissipation for several vibrational multipolarities. For the considered 
oscillation modes an increasing divergence of the nucleonic trajectories from the adiabatic to the 
resonance regime was observed. Also, a peculiar case of intermittency is reached in the vicinity of the 
resonance, for the monopole case. We examine the order-to-chaos transition by performing several 
types of qualitative analysis including sensitive dependence on the initial conditions, single-particle 
phase space maps, fractal dimensions of Poincare maps and autocorrelation functions. 

PACS numbers: 24.60.Lz, 05.45.-a, 05.45.Pq, 21.10.Re 



> 
o 

O 



I. INTRODUCTION 

Deterministic chaos, is usually defined as irregular, un- 
predictable behaviour of the trajectories generated by 
nonlinear systems whose dynamical laws, involving no 
randomness or probabilities, predict a unique time evo- 
lution of a given system. 

Over the last two decades an increasing number of pa- 
pers have treated the study of the deterministic chaotical 
behaviour of Fermi nuclear systems pl-(27l| . The interest 
for analyzing the order-to-chaos transitions on such sys- 
tems was linked to the problem of the onset of dissipation 
of collective systems through mainly one-body and two- 
body processes. Among these we mention the interaction 
of the nucleons with the potential well, the evaporation 
of individual nucleons in nuclear peripheral interactions, 
and the collisions between nucleons without taking into 
account the Pauli blocking effect. 

These kinds of analyses were performed for the first 
time by Burgio, Baldo et al. considering a system 

of nucleons which move within a container modelled as a 
Woods-Saxon type potential and kick the container walls 
with a specific frequency. They discuss the damping of 
the movement and the relation with order-chaos transi- 
tion in single-particle dynamics. 

On that classical model worked Papachristou and col- 
lab. [2^, studying the decay width of the Isoscalar Giant 
Monopole Resonance for various spherical nuclei. Fol- 
lowing also that formalism, the beginning of the chaotic 
behaviour for a number of nucleons in various dynamical 
regions at several multipolarities was surveyed [29l433l | . 



In this paper, we investigated the chaotic behaviour of 
a single-nucleon in a two-dimensional {2D) deep Woods- 
Saxon potential well for specific physical phases. A qual- 
itative picture of the achievement of deterministic chaos 
was shown for a comparative study between the adiabatic 
and the resonance stage of the nuclear interaction. 

Close to resonance we obtained characteristics of the 
intermittency regime, i.e. sudden change to a laminar 
behaviour (so-called intermission) of a specific signal be- 
tween two turbulent phases, which has been detected over 
and over again in a plethora of experiments regarding the 
Rayleigh-Benard convection, the driven nonlinear semi- 
conductor oscillator, the Belousov-Zhabotinskii chemical 
reaction, and the Josephson junctions (for e.g., [ssl - li^ ). 

Albeit intermittency is a well-known phenomenon for 
billiards [isl - ldH , in particular for Hamiltonian ^stems 
with divided phase space (e.g., mushroom 49-53] and 
annular billiards [53|). and for connected Hamiltonian 
systems [s^ , we showed that this property also holds for 
a Woods-Saxon billiard container with inelastic particle- 
wall interactions. 



II. BASIC FORMALISM 

We chose for the present analysis a simple dynamical 
system as considered in several previous papers by Burgio 
and coUab. [J-Q- This system contains a number of A 
spinless and chargeless nucleons, with no internal struc- 
ture. The nucleons move in a 2D deep Woods-Saxon 
potential well considered as a "nuclear billiard" and hit 
periodically the oscillating surface of the well with a cer- 
tain frequency. The Bohr Hamiltonian of such a system 
in polar coordinates is considered as: 
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where {prj,Pej,Pa} are the conjugate momenta of the 
particle and coUective coordinates {vj, 9j, a}, the nucleon 
mass m is 938 MeV, ft is the oscillating frequency of the 
collective variable a, and M = tuARq is the Inglis mass. 

The Woods-Saxon potential is constant inside the bil- 
liard and a very steeply rising function on the surface: 



V(r,,R{e,)) = 



Vo 



1 + exp 



rj-R(9j,a) 



(2) 



with Vb = —1500 MeV, deep enough to prevent the 
escape of the nucleons regarded as classical objects for 
the present analysis. For the same reason, the diffusivity 
coefficient a has a very small value 0.01 fm. The vibrat- 
ing surface can be written as in [J-Q, depending on the 
collective variable and Legendre polynomials Pl (cos 6j ) : 



R,=R( 



= i?o [1 + oiPl {cos0j] 



(3) 



where i?o = 6 fm, and L the multipolarity vibration 
degree of the potential well is for the monopole, 1 for 
the dipole, and 2 for the quadrupole case. 

Once the hamiltonian is chosen, the numerical simula- 
tions are based on the solution of the Hamilton equations: 
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The Hamilton equations were solved with a Runge- 
Kutta type algorithm (order 2-3) having an optimized 
step size and taking into account that the absolute error 
for each variable is less than 10~^. Total energy was 
verified to be conserved with high accuracy to a relative 
error level of 10~*. 

The equilibrium deformation parameter a, which is the 
mean collective variable, can be calculated (for e.g., for 
L = 0) by equating the mechanical pressure of the wall, 

Pwall • 
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The pressure exerted by the particles is Ppart, P de- 
notes the particle density, and T is the apparent temper- 
ature of the system that equals the 2D kinetic energy, 
using the natural system of units {h ^ c = ks = 1): 



Ppart pT-^ 
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Thus, one gets the equation for the equilibrium value 
of the collective coordinate in 2D, in the monopole case: 
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2T 
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(9) 



Then a small perturbation of this collective variable 

was considered a{t = fm/c) = a -1-0.15 [iHll and the 
evolution of the physical system was thoroughly investi- 
gated. 



III. THE QUALITATIVE ANALYSIS OF THE 
ROUTE TO CHAOS 



A. On the resonance condition 

One can choose the wall oscillation taking place close to 
adiabatic conditions, imposing a wall frequency smaller 
than the single-particle one. Thus, the frequency of vi- 
bration flad was chosen less than 0.05 c/fm, which corre- 
sponds to an oscillation period equal to: 



27r 

Twaii = 7r—> 125.66 fm/c. 
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By introducing the maximum particle speed: 
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and the parameters as in [H-i]: i?o = 6 fm and T = 
36 MeV, one can obtain the value for the single-particle 
period: 



Tpart = — ^ ~ 43.33 fm/c. 

V 



(12) 



In addition to [ll-Q we introduced a physical constraint 
to this elementary physical system and continued that 
type of analysis necessary for the study of a nonintegrable 
dynamical system. At the beginning we considered a 
physical situation and we chose instead a static vibrating 
"nuclear billiard", a projectile nucleus having the same 
properties, colliding with a target nucleus. It is well- 
known that the nuclear interaction, at incident energies 
ranging from MeV to GeV, can result in a multitude 
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of processes from the nuclear evaporation to complete 
fragmentation or multifragmentation, according to the 
impact parameter. 

It was shown in [sgI . [57| that during this kind of pro- 
cesses even for peripheral events an unnegiigeable amount 
of energy is transferred by nucleon-nucleon scattering to 
the nucleons of the projectile and not only the trans- 
verse momentum distributions, but also the longitudinal 
momentum distributions as measured in the projectile 
fragmentation rest frame can reveal the centrality status 
of the interaction. It can also offer a hint on the apparent 
temperature of a Fermi gas of nucleons which was found 
to be 1571 n ear the isotopic temperatures, i.e. several 
MeV [HileH. 

It was therefore supposed that the target fragmenta- 
tion can be associated with a resonance process. In or- 
der to obtain such behaviour, the wall frequency was 
gradually increased to the resonance frequency fires = 
0.145 c/fm. However, nuclear evaporation or plain 
breakup of a projectile nucleus can take place long before 
this regime is achieved by redistributing energy between 
the nucleons themselves and also between single-particle 
degrees of freedom and collective ones. Individual nucle- 
ons or clusters can thus have enough kinetic energy to 
escalade the wall barrier. 

We should also emphasize that we can either have the 
case that can be put in correspondence with a nuclear 
collision process, i.e. the variation of the nuclconic fre- 
quency oscillation as the apparent temperature of the 
nucleons in the nuclei increases (Eqs. 11 and 12), main- 
taining the potential well vibration constant, or respec- 
tively, the inverse situation in which the period between 
two consecutive collisions of the nucleon with the self- 
consistent mean field is kept invariable, while modifying 
the oscillation modes of the nuclear surface. The latter 
regards our studied case and is the reversed physical case 
previously described. It was used because of the specific 
choice of the "toy model" parameters described in [H-Q. 

But the most realistic evolution of the nucleons in a 
chosen potential can assume a simultaneous variation of 
both angular frequencies. The resonance condition of 
the coupled classical oscillators should remain however 
an important condition for a rapid appearance of a de- 
terministic chaotical behaviour of the physical system in 
study at different time scales. A proper analysis of a 
system should provide the variation of the collision ra- 
dian frequency of the nucleons inside the "billiard" as 
the apparent temperature increases and the change in 
the vibrating potential period, supposing that the multi- 
polarity increases when pumping energy in the "nuclear 
reservoir" during interaction. We can for example use in 
simulations, for nuclei with a large number of nucleons, 
the Liquid Drop Model or the Collective Model, which 
predict a frequency of vibration as function of the multi- 
polarity deformation degree: 



with Cl being the elasticity coefficient, and Bl the 
mass coefficient for the oscillator of L multipolarity. 

We will briefly discuss on how the resonance condition 
might look like and for which particular case(s) it can be 
applied. By reverting to the set of differential nonlinear 
Hamilton equations and combining the last two of them 
(see Eq. 6), it emerged: 
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One can now compare the resulted equation with a 
harmonic oscillator system in the presence of dissipative 
processes ("damped oscillator") and of a external har- 
monic driving force of _Fb amplitude: 



-r a +n'^a = 



FoRo 
M 



■ cos Lot. 



(15) 



In a first approximation we can neglect the damping 
constant F, as noticed from the Figs. [T]|4l Besides, dis- 
sipation will appear if only the collective coordinate a 
is averaged over a large number of events, for all mul- 
tipolaritics considered P, [1]. For the uncoupled equa- 
tions (UCE) case there is no damping either, as the so- 
lution of a homogeneous Eq. 14 is a purely harmonic one: 
a — ao cos fit. 

Also, we remark that for the dipole collective oscil- 
lations: Pi (cosOj) = cos LUjt and that, for a statis- 
tical ensemble of A nucleons, the one-particle radian 
frequency can be approximated with a constant, ujj sa 
uj. For the Woods-Saxon potential energy: V{t) = 
^f=iV{rj,6j,a), it can always be derived a conserva- 
tive force Fq: 



Foit) = -J2 



dr-i 



j=i 1 



(16) 



(13) 



A similarity between Eq. 14 and 15 was eventually 
obtained, because once shifting towards collectives of nu- 
cleons, the total force becomes a sum over a number of 
time dependent forces, each acting on a single nucleon. 
Therefore, the net force for a pack of nucleons tends to 
vary little with time: Fq ^ FQ{t), and the tendency in- 
creases with A. 

In conclusion, the resonance condition, as stated in Eq. 
14, can be definitely used for the multi-particle couplings 
with a dipolar shape of the potential well. This was ver- 
ified by studying the variance of the largest Lyapunov 
exponent and of the Kolmogorov-Sinai entropy with the 
radian frequency of wall oscillation on a ten- nucleons sys- 
tem [3^. For the other multipole degrees taken into ac- 
count, and as well, for the dipole oscillations coupled 
with a single-nucleon dynamics, the nonlinear character 
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of the differential Equation 14 offers a somewhat intri- 
cate perspective over obtaining an analytical generalized 
resonance condition, and for the moment is left aside for 
further investigation. 

As we confine the following analysis on the single- 
particle chaotical dynamics, the frequency matching of 
the coupled oscillators: uipart = ^res = 0.145 c/fm, will 
be generically designated as the resonance stage of the 
interaction. 



B. Dependence on the initial conditions 

In order to detect chaos in simple systems, several cus- 
tomary methods are used (6^ . The sensitive dependence 
on the initial conditions, the so-called "Butterfly effect" 
is the first type of analysis presented in this article. For 
a given multipolarity degree we studied the time depen- 
dence of both, single-particle and collective variables, 
with small perturbations applied to a single parameter 
(for e.g., Ar = 0.01 fm). This type of analysis is pre- 
sented for four cases, ranging from UCE to quadrupole 
deformations of the two-dimensional wall surface, and for 
the two physical regimes in study, adiabatic and of reso- 
nance, respectively (Figs. [l]|4]). 

At a first glance, the time dependence of the col- 
lective degrees of freedom "looks chaotic" for all con- 
sidered situations with the exception of (4 -I- 2) uncou- 
pled nonlinear differential equations case (Fig. [T]), as 
one would naturally expect. In this case, a periodical 
structure would definitely appear in sharp contrast with 
the quasi-periodicity shown for different multipole col- 
lective modes. Also, the time evolution of the single- 
particle variables is clearly chaotic in all cases where the 
strong coupling between the (4 -I- 2) nonlinear differential 
equations describing the single-particle dynamics appears 

(Figs. mM- 

The sensitive dependence on the initial conditions, i.e. 
the decoupling of the trajectories at macroscopic scale, is 
found to give the first hints on the behaviour of the nu- 
clear system, in its evolution from the quasi-stable states, 
which can hardly develop a chaotic motion in time (adi- 
abatic state), to the unstable ones, characterized by a 
rapid divergence of the particle trajectories in the phase 
space (resonance regime). Also, in the monopole case 
the order-to-chaos transition clearly shows an intermit- 
tent pattern at = 0.1 c/fm (Fig. [5]). 

C. Maps of phase space 

Another type of qualitative analysis, indicating differ- 
ent ways toward a chaotic behaviour is based on the one- 
dimensional maps of phase space. 

Thus, for a temporal scale of 3, 200 fm/c we repre- 
sented the evolution of the multipolarity deformation de- 
gree of the potential well in the phase space of the vi- 
brational variables (a <-^- Pa), as well as in the ID phase 



space of a single-particle (r <->Pr)- This type of analy- 
sis was performed for the UCE, monopole, dipole and 
quadrupole cases, using the same physical conditions, 
from the adiabatic to the resonance phase (Figs. [5l fT^ . 

A remark that can be formulated from the study of 
these phase space maps is related to the form of the tra- 
jectories specific to those described by a harmonic oscil- 
lator, especially for the (a <-> Pa) maps. It can be easily 
put to the test (Figs. El [71 [3 and [Tl|) that the ellipse 
area Si-„ is proportionally inverse with the wall vibration 
frequency: 

27r ■ Ecou , . 

where Ecoii is the sum of the last two terms from (1), 
being the energy of the collective nucleonic motion. 

The phase space filling degree raises as the vibration 
wall frequency moves towards the resonance value for all 
multipolarities took into account (Figs. [Tl fT^ . with the 
exception of the uncoupled nonlinear differential equa- 
tions case (Figs. [5] and [6]). This confirms that the char- 
acteristic time of the macroscopic decoupling of the nu- 
cleon trajectories in the Woods-Saxon potential evolves 
towards smaller values, once D, is increased and also 
that the coupling between the collective variable motion 
and the particle dynamics is essential in amplifying the 
chaotic behaviour. 

An intriguing aspect is revealed by the comparison of 
the filling degrees of the phase space as a function of 
the nucleon oscillation frequency in the chosen potential. 
One would expect that the trajectories degenerate from a 
simple orbit towards a compact filling of the (a <-> Pa) or 
(r <-> Pr) plane. Moreover, we observe several attraction 
basins around a few standard orbits, preeminently in the 
monopole case at = 0.1 c/fm, characteristic feature of 
an intermittent behaviour (Figs. |7land|8]). 

D. Fractal dimensions of Poincare maps 

The one-dimensional phase space maps offer a hint on 
their filling degree in time. In the same manner one can 
use the Poincare maps for the " nuclear billiard" regarded 
as a deterministic dynamical system. When analyzing 
the behaviour of close trajectories in the phase space, 
starting from a periodical orbit, Poincare showed [63i, ,6^ 
that there are only three distinct possibilities: the curve 
can be either a closed one (with a fix distance to a pe- 
riodical orbit), or it can be a spiral that asymptotically 
wraps/unfolds to/from the periodical trajectory. 

In order to simplify the temporal analysis for peculiar 
orbits, Poincare suggested a simple and effective method. 
By choosing a transversal section with (TV — 1) D that 
intersects a ND geometrical variety, the sequence of the 
intersection points is easier to be studied than the whole 
ND curve. Therefore, the cases previously described for 
the spiral generate a series of points that draw near in or 



r £ [2.;2.04] fm; A r = 0.01 fm; N , = 1 ; Uncoupled Eqs. 



nucleons 

n = 0.010 c/1m n _ = 0.050 c/lm fl = 0.100c/1m 



12 = 0.145 c/1m 






200 400 600 800 200 400 600 800 200 400 600 

200 ....... 200 r, 200, 




200 400 



200 400 



MMA ^ MMA ^ V\m/v ■ \Am/v 



262 
2G1 
2G0 
269 



200 400 600 800 200 400 600 800 200 400 600 

262 

M- 

260 ■ 



200 400 



259 



262 
261 
260 



259 



200 400 600 800 200 400 600 800 200 400 600 



262 

2G0 
259 



200 400 






200 400 600 800 200 400 600 800 200 400 600 



200 400 600 800 





200 400 600 800 200 400 600 800 200 400 600 

1 (fm/c) 1 (1m/c) t (fm/c) 



200 400 600 800 

t (1m/c) 



FIG. 1. The sensitive dependence on the initial small perturbation of the radius parameter (0.01 fm/c) when adiabatic and 
resonance conditions are imposed (uncoupled one-nucleon and collective degrees of freedom case). 



move away from the intersection point of the periodical 
orbit with the transversal section. 

Applying the Poincare maps type analysis to the sim- 
ple chosen physical system, one would expect that the 
compact filling up of a IZ? map of phase space for vibra- 
tion frequencies of the Woods-Saxon well that match the 
one-particle oscillation frequencies inside the "billiard" 
(Figs. [5l fT2|) . to reverberate by a large density of chaotic 
distributed intersection points. 

In order to estimate the degree of fractality of such 
densities, we computed the fractal dimensions d/ of the 
Poincare maps with one-nucleon radial degrees of free- 
dom (r <-> Pr)) when choosing for the transversal section 
the polar pair {9 pe) as a constant. As only for the 
monopole and UCE cases the orbital kinetic momentum 
is a constant of motion, we plotted the one-particle ra- 
dius as a function of radial momentum when the polar 
coordinate is not kept constant, but quasi-constant in or- 
der to increase the probability of intersecting the section 
for a given time evolution of the system: 



(18) 



We then chose a very small value for the (6* ^ pe) 
thickness: eg = 10~^ radians and let the system evolve 
for a period of At = 5 ■ 10® fm/c. 

The values were calculated using the box counting al- 
gorithm in a specific Visual Basic 6 application [gl, 
on Poincare maps. The fractal dimension of a map cov- 
ered by N boxes of length r is described by the following 
expression (see, for e.g., (67l[68j): 



log N (r) 
r^b log (1/r) 



lim 



(19) 



where we took: r = 1/2", < n < Umax- The maxi- 
mum value Umax denotes the highest chosen resolution. 

N(r) 

We calculated d/ as slopes of the linear fits of logj 

1 /r 

versus n = logj points, where we took the inferior limit 
of log{l/r) to be 5, and the superior threshold according 
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FIG. 2. The sensitive dependence on the initial small perturbation of the radius parameter (0.01 fm/c) when adiabatic and 
resonance conditions are imposed (monopole case). 



to the highest resolutions for which a saturation plateau 
appears: N N{r) (Fig. [T5|) . At resolutions above 
threshold, no additional information on the points distri- 
bution can be gained. 

The points from the Poincare maps arc distributed in 
such a way that up to a certain resolution (32 x 32 pixels) 
the algorithm does not take into account the individual 
pixels. In this first region analyzed the alignment is quite 
good, denoting a correlation between points. The infor- 
mation is not uniformly distributed in the phase space 
but represents a "cloud" of points with df greater than 
1. This large-scale region thus offers a global image of 
the Poincare maps. 

Above 32 x 32 resolution appear the effects linked with 
the contribution of the individual dots of the Poincare 
maps, reflected by significant variations of the fractal di- 
mensions. By choosing the points for the fitting proce- 
dure in such a way, we are also in agreement with the 
following criterion: the fitted points should be singled 
out so their associated resolutions are to be pertained 
to at least three successive decimal logarithmic intervals 



TABLE I. The d/ of the (r 
maps 



Pr) Poincare single-particle 



Oscillation frequency Uncoupled eqs. Monopole case 



J7 ad = 0.010 c/fm 0.1194 

f2 ad = 0.050 c/fm 0.1666 

Q =0.100 c/fm 0.1840 

f^res = 0.145 c/fm 0.2227 



0.1262 
0.1870 
0.1790 
0.2010 



[69l |. In Figure [13] we selected two points with resolutions 
belonging to the first interval [lO 10^] , three points to 
[10^ ^ 10^] , and we took one point in the [lO^ lO"] 
range, just until constant values N are reached. 

The small-scale df thus computed (Table |T| represent 
the filling degrees of a fine detailed phase space and can 
be correlated with the Shannon entropies calculated in 

It can be remarked that the fractal dimensions gener- 
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FIG. 3. The sensitive dependence on the initial small perturbation of the radius parameter (0.01 fm/c) when adiabatic and 
resonance conditions are imposed (dipole case). 



ally increase with the specific vibration frequencies, from 
the quasi-stationary equilibrium regime at adiabatic os- 
cillations to the unstable chaotical states of the nucle- 
onic system in the resonance domain. At the same time, 
it should be mentioned that for the monopole case the 
fractal dimension has an intermittent deportment close 
to the resonance phase of interaction, ai Q = 0.1 c/fm. 

E. Autocorrelation functions 

The transition toward a quasi-periodical, aperiodical 
or chaotic behaviour can be also analyzed with the au- 
tocorrelation function of the variables characteristic for 
the temporal evolution of the studied physical system. 
When the "nuclear billiard" evolves to chaos, essential 
changes can be noticed in the shape of the autocorrela- 
tion function of a specific variable (continuous or discrete 
distributed). This type of function measures the correla- 
tion between a sequence of signals, and is usually defined 
as: 



C(t)= hm - / Ht)-Ht + r) dt; (20) 



1 f 

= x{t)- lim - / x{t) dt; (21) 



C(r)= hm -^^(i„).e(i„+i); (22) 



? {tn) = X (i„) - lim - ^ i (t„) . (23) 

For a laminar regular regime this function either has a 
constant value or presents decreasing oscillations in time. 
During a chaotic phase, it falls in, having an exponen- 
tially decreasing behaviour for uncorrelated signals, as 
shown in [Tlj . 
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FIG. 4. The sensitive dependence on the initial small perturbation of the radius parameter (0.01 fm/c) when adiabatic and 
resonance conditions are imposed (quadrupole case). 



We represented the autocorrelation function for the 
study of the single-particle (Fig. and collective dy- 
namics (Fig. [T5)) . The analysis was thus done auto- 
correlating the radial momentum p,. of the particle and 
respectively the collective variable a for the UCE case 
and for all multipolarities, from all stages of interaction 
taken into account. One conclusion, that can be easily 
drawn from these figures, confirmed the previous results 
obtained with the other two methods. Thus, it can be 
noticed the increase of the chaotic bearing once the os- 
cillation frequency of the potential well is varied from 
the adiabatic to the resonance regime, for all degrees of 
multipolc. 

When studying the case with single-particle degrees of 
freedom uncoupled from the collective ones, the domi- 
nant behaviour is aperiodical oscillation in time charac- 
teristic for a steady adiabatic phase of the interaction, 
even though the vibration frequency is varied. Only at 
resonance the shape is somewhat changed, indicating a 
greater degree of chaos. 

As for the monopolc to quadrupole deformations the 



decreasing of the autocorrelation function is indeed of 
the exponential type, steeper as the radian frequency of 
oscillation is raised to fires = 0.145 c/fm. 

The exception is again found in the monopole case at 
= 0.1 c/fm frequency of oscillation, when the inter- 
mittent phase of the interaction, reflected by aperiodical 
oscillations, points out a steady behaviour prior to the 
resonance regime. 

IV. CONCLUSIONS 

A comparative study was done between the interest- 
ing physical regimes of nuclear interaction: adiabatic and 
resonance, giving at this level only a qualitative picture 
of the possible scenarios towards a pure deterministic be- 
haviour of chaotic type of the studied nucleonic system. 

We envisaged the single-nucleon dynamics in a Woods- 
Saxon potential. The coupling between individual and 
collective degrees of freedom was shown to generate dif- 
ferent paths to chaos, according to the order of multi- 
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{a - p ) phase-space map; N = i ; a t = 3200 fm/c; Uncoupled Eqs. 
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;a - p ) phase-space map; N = 1 ; A t = 3200 fm/c; L = (Monopole) 
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FIG. 5. The phase space of the collective degrees of freedom FIG. 7. The phase space of the collective degrees of freedom 
{a — Pa) for different wall frequencies (uncoupled single and (a — Pa) for different wall frequencies (L — 0). 
collective degrees of freedom case) . 
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FIG. 6. The phase space of the single-particle degrees of free- 
dom {r — pr) for different wall frequencies (uncoupled single 
and collective degrees of freedom case) . 



polarity. In general, for all vibrational modes at the 
resonance frequency of oscillation, the onset of chaoti- 
cal behaviour was found to be earlier than at any other 
adiabatic vibrations of the 2D potential well. 

Also, a phase alternation of periodical and chaotic dy- 
namics was found in the monopole case of nuclear wall 
oscillation at 17 = 0.1 c/fni, revealing a laminar dynamics 
prior to the resonance stage of interaction. 

In order to verify the aforementioned results, the study 
was completed with the inclusion of several quantitative 
analyses, inter alia we mention: power spectra. Shannon 
informational entropies, and Lyapunov exponents [7Q] . 



FIG. 8. The phase space of the single-particle degrees of free- 
dom (r — Pr) for different wall frequencies (L = 0). 
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FIG. 9. The phase space of the collective degrees of freedom 
{a — Pa) for different wall frequencies (L — 1). 
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(r - p ) phase-space map: N = i ; a t = 3200 fm/c: L = 1 (Dipole) 



Uncoupled Eqs. 
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FIG. 10. The phase space of the single-particle degrees of 
freedom {r — p,-) for different wall frequencies (L = 1). 



(a - p ) phase-space map; N . = 1 ; A t = 3200 fm/c; L = 2 (Quadrupole) 
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FIG. 13. The fractal dimensions of the (radius ^ radial mo- 
mentum) Poincare one-nucleon maps as linear slopes. 
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FIG. 11. The phase space of the collective degrees of freedom 
(a — pa) for different wall frequencies (L = 2). 
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FIG. 12. The phase space of the single-particle degrees of 
freedom (r — Pr) for different wall frequencies (L = 2). 



ii SIS 
I I, (I 



T (fm/c) 



T (fm/c) 



T (fm/c) 



T (Im/c) 



FIG. 14. The autocorrelation function which correlates the 
radial momentum of a single-particle in time for various stages 
of nuclear interactions (from adiabatic to resonance regimes) 
and for several nuclear multipole deformations. 
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